setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller")
m = read_csv("Mueller.csv")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller")
m = read_csv("Mueller.csv")
library(ineq)
library(readr)
library(doSNOW)
library(foreach)
library(snow)
library(tcltk)
library(data.table)
library(extremevalues)
library(truncnorm)
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller")
m = read_csv("Mueller.csv")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller")
m = read_csv("Mueller.csv")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Level Inequality Fit 2")
p = read_csv("Mueller Mean and Gini by Level.csv")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Level Inequality Fit 2")
m = read_csv("Mueller Mean and Gini by Level.csv")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Level Inequality Fit 2")
m = read_csv("Mueller Mean and Gini by Level.csv")
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Level Inequality Fit 2")
d = read_csv("Mueller Mean and Gini by Level.csv")
sigma = 2*invErf(d$gini)
setwd("/home/blair/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Level Inequality Fit 2")
d = read_csv("Mueller Mean and Gini by Level.csv")
# actual data
###########################################
sigma = 2*invErf(d$gini)
sigma = 2*invErf(d$gini)
mu =  log(d$Mean) - 0.5*sigma^2
n = round(d$N*10)
m = data.table(sigma, mu, n)
m = na.omit(m)
sigma = rep(m$sigma, m$n)
mu = rep(m$mu, m$n)
n = length(mu)
x = rlnorm(n, mu, sigma)
ineq(x)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Test Economy/Test Economy.R')
ineq(x)
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Test Economy/Test Economy.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Test Economy/Test Economy.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Test Economy/Test Economy.R')
source('~/Desktop/Empirical Research/Income Distribution/Gini Ratio/Hierarchical Levels/Empirical/Aggregate Studies/Mueller/Test Economy/Test Economy.R')
